home *** CD-ROM | disk | FTP | other *** search
/ NASA Climatology Interdisciplinary Data Collection / NASA Climatology Interdisciplinary Data Collection - Disc 4.iso / software / grads / lib / clhilo.gs < prev    next >
Text File  |  1998-04-23  |  2KB  |  71 lines
  1. function clhilo(args)
  2. *
  3. *    plot H/L using Steve Lords' algorithm
  4. *
  5.   var=subwrd(args,1)
  6.   maxmin=subwrd(args,2)
  7.   fmt=subwrd(args,3)
  8.   cint=subwrd(args,4)
  9.   rad=subwrd(args,5)
  10.  
  11.   if(maxmin = "maxmin" | maxmin = "MAXMIN") ; mm=0.0 ; endif
  12.   if(maxmin = "max" | maxmin = "MAX") ; mm=1.0 ; endif
  13.   if(maxmin = "min" | maxmin = "MIN") ; mm=-1.0 ; endif
  14.  
  15.   if(fmt='') ; fmt=i5; endif
  16.   if(rad = ''); rad=500.0; endif
  17.  
  18.   if(cint = '-999' | cint = '') 
  19.     'set gxout stat'
  20.     'd 'var
  21.     dum=sublin(result,9)
  22.     say "dum="dum
  23.     cint=subwrd(dum,7)
  24.     say "cint set to "cint
  25.     'set gxout contour'
  26.   endif
  27.  
  28.   udfile='udf.clhilo.vals'
  29.   '!rm 'udfile
  30.   say "rad,cint="rad","cint
  31. * 
  32. *    call the udf to calc H/L
  33. *
  34.  'd clhilo('var','mm','rad','cint','fmt')'
  35. *
  36. *    plot the H/L
  37. *
  38.   while (1)
  39.     res=read(udfile)
  40.     rc=sublin(res,1)
  41.     if (rc != 0)
  42.       if(rc = 1); say "open error for file "udfile; return; endif
  43.       if(rc = 2); say "end of file for file "udfile; break; endif
  44.       if(rc = 9); say "I/O error for file "udfile; return; endif
  45.     endif
  46.     dum=sublin(res,2)
  47.     maxormin=subwrd(dum,1)
  48.     lon=subwrd(dum,2)
  49.     lat=subwrd(dum,3)
  50.     val=subwrd(dum,4)
  51.  
  52.     'q w2xy 'lon' 'lat
  53.     x=subwrd(result,3)
  54.     y=subwrd(result,6)
  55.     if(maxormin = "max"); str="H"; 'set string 4 c 6'; endif
  56.     if(maxormin = "min"); str="L"; 'set string 2 c 6'; endif
  57.     'set font 5'
  58.     'set strsiz 0.1'
  59.     'draw string 'x' 'y' 'str
  60.     'set font 0'
  61.     'set strsiz 0.05'
  62.     'set string 1 tc 6'
  63.     ytop=y-0.1
  64.     'draw string 'x' 'ytop' 'val
  65.     
  66.   endwhile
  67.   rc=close(udfile)
  68.   say "End of clhilo"
  69.  
  70. return
  71.